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ABSTRACT 

We study the low-contrast Fell emission blends in the ultraviolet (1250- 
2200 A) and optical (4000-6000 A) spectra of the Seyfert 1 galaxy NGC 5548 
and show that these features vary in flux and that these variations are correlated 
with those of the optical continuum. The amplitude of variability of the optical 
Fell emission is 50% - 75% that of H/5 and the ultraviolet Fell emission varies 
with an even larger amplitude than H/5. However, accurate measurement of the 
flux in these blends proves to be very difficult even using excellent Fe ll templates 
to fit the spectra. We are able to constrain only weakly the optical Fe ll emission- 
line response timescale to a value less than several weeks; this upper limit exceeds 
all the reliably measured emission-line lags in this source so it is not particularly 
meaningful. Nevertheless, the fact that the optical Fell and continuum fiux vari- 
ations are correlated indicates that line fiuorescence in a photoionized plasma, 
rather than collisional excitation, is responsible for the Fe ll emission. The iron 
emission templates are available upon request. 

Subject headings: galaxies: active - galaxies: nuclei — galaxies: Seyfert — 
quasars: emission lines 

1. INTRODUCTION 

Emission from singly ionized iron is one of the principal coolants of the broad-line region 
(BLR) in active galactic nuclei (AGNs): indeed, the prominent blend of Fell and Balmer 
continuum emission in the near UV spectrum, a feature often referred to as the "small blue 
bump," produces as much as one third of the line fiux in some objects (Wills, Netzer, & 
Wills 1985; Maoz et al. 1993). Remarkably, however, the characteristics of Fell emission in 
AGNs are very poorly understood from both theoretical and observational perspectives: 
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1. The excitation mechanism that produces Fell emission in AGN spectra is not known. 
Theoretical models, based on either photoionization or collisional excitation, have for 
a long time fallen short in reproducing the observed Fell emission (e.g., Collin & Joly 
2000; Baldwin et al. 2004). Also, conflicting results appear in the literature as to 
the importance of the hard X-rays as an excitation agent as originally anticipated 
by Davidson & Netzer (1979) and as to that of the slope of the soft X-ray spectrum 
(e.g., Wilkes, Elvis, & M^Hardy 1987; Shastri et al. 1994; Boroson 1989; Zheng & 
O'Brien 1990; see also Wilkes et al. 1999). Even with much-improved current models, 
either significant microturbulence (or velocity shear) in a photoionized medium or a 
contribution from a coUisionally ionized gas seems to be required to account for the 
strength and shape of Fell emission in AGN spectra (e.g., Baldwin ct al. 2004). The 
conditions that are necessary to reproduce the observed Fell spectra do not provide 
clear and straightforward diagnostics to distinguish between these possibilities. 

2. Strong, broad Fell emission pervades the UV through IR spectra of AGNs, contami- 
nating other spectral features. The complex atomic structure of Fe"*" results in tens of 
thousands of transitions (e.g., Verner et al. 1999; Sigut & Pradhan 2003) that produce 
blended features across the UV, optical, and IR spectrum. Isolating individual emis- 
sion lines under these circumstances is challenging, leading to systematic uncertainties 
in measurement of emission-line fluxes, widths, and proflles. The intrinsically large 
line widths of AGN emission lines add to the complexity of this issue. 

Further underscoring the importance of understanding AGN Fe ll emission is the possibility of 
using AGN spectra to measure the cosmic evolution of the abundance of iron, thus providing 
potentially strong constraints on the history of chemical enrichment due to star formation 
(e.g., Hamann & Ferland 1999; Dietrich et al. 2003 and references therein, but see also 
Verner et al. 2004; Baldwin et al. 2004). 

A potentially powerful way of exploring the Fe ll emitting region in AGNs is through flux 
variabihty, specifically by emission-line reverberation mapping (Blandford & McKee 1982; 
Peterson 1993, 2001). At minimum, reverberation mapping can be expected to yield the 
typical distance of the Fe ll-emitting gas from the central continuum source through mea- 
surement of the time-delayed response of the emission lines to continuum variations. This 
would allow us to constrain observationally some of the physical conditions conducive to 
strong Fell emission. For example, variability time scales for the optical and UV Fell emis- 
sion, respectively, will help discern whether or not the multiplets at these energies originate 
in a common region. This is a key issue for understanding their relative intensity ratio in 
terms of current Fell models which fail to reproduce observations (Baldwin et al. 2004). At 
the present time, studies of Fe ll fiux variability have been very limited: 
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1. Maoz et al. (1993) used UV data obtained with the International Ultraviolet Explorer 
(lUE) and optical data obtained with multiple ground-based telescopes by the Inter- 
national AGN Watch^ over an 8-month period in 1989 to study the variability of the 
small blue bump in NGC 5548. They concluded that the small blue bump emission 
arises approximately at the same distance as the Lya emission, about 10 light days 
from the continuum source. While in many ways, this represents the best study to 
date on variability of the small blue bump, it must be noted that (a) the quality of 
lUE data in this part of the spectrum (i.e., obtained with the LWP camera) is quite 
poor and (b) there are relatively few epochs in which there is no gap in wavelength 
coverage between the UV and optical spectra, making measurements of the small blue 
bump variability systematically somewhat uncertain. 

2. Goad et al. (1999) observed the small blue bump in NGC 3516 on six occasions over 
an 11-month period in 1995-96. Remarkably, while the continuum and high-ionization 
emission hnes varied in flux by a factor of 2, the Mgll A2800 and Fell UV emission did 
not vary by more than 7% over this time. 

3. KoUatschny, Bischoff, & Dietrich (2000) carried out a 20-year study of emission-line 
variations in NGG 7603 and found that the optical Fe ll blends vary with an amphtude 
comparable to that of the Balmer hnes. 

4. KoUatschny & Welsh (2001) searched for and found no optical Fell variability in Mrk 
110 to within their sensitivity of 10%, although clear optical Fell variability is seen in 
the case of Fairall 9 with an amplitude half that of H/? (KoUatschny & Fricke 1985). 

5. Sergccv ct al. (1997) attempt to constrain the time delay of the Fell multiplets in the 
wings of [Olll] A5007 in the spectrum of NGC 5548. They find a time delay of a few 
hundred days. There arc, however, a number of serious difficulties with this result 
arising from strong blending of lines, data quality, and time-series problems. 

6. A few other studies make claims of Fe II variability (KoUatschny ct al. 1981; Salamanca, 
Alloin, & Pelat 1995; Doroshenko et al. 1999). On the basis of spectra of narrow-line 
Seyfert 1 (NLSl) galaxies taken a year apart, Giannuzzo & Stirpc (1996) find that the 
optical Fell does seem to vary in these sources. For the three sources with detected 
variability in both H/3 and Fell, the Fell/H/3 flux ratio ranged from 0.25 to 1.0; this 
would put Fe II variability below the detection threshold for the sources in which it was 
not detected. 
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Given the importance of understanding Fell emission in AGNs and given the hmited 
attention this topic has received in the past, we have undertaken an examination of Fell 
variabihty in NGC 5548. In terms of emission-hne variabihty, NGC 5548 is by far the best- 
studied AGN (see Peterson et al. 2002 and references therein). On the other hand, as Baldwin 
et al. (2004) point out. Fell emission is not especially strong in this source and other sources 
may be better suited to such a study. While we agree with this statement, we also feel that 
the wealth of data on other lines in NGG 5548 affords a special context in which to examine 
Fell variabihty and, of course, suitable data for this study already exist. Moreover, cursory 
consideration of Fell variability in NGG 5548 yields mixed impressions. For example, the 
difference between spectra taken at recent historical maximum and minimum flux states (see 
Peterson ct al. 1999, Figure 3) seems to clearly show optical Fell emission, indicating that 
these blends do indeed vary in this source. However, Kollatschny & Welsh (2001) attribute 
the apparent variations solely to changes in HeIlA4686. If indeed the optical Fell emission 
varies with time, then we ought to see evidence of these features in the difference or rms 
spectra (defined below in §2.3.2.) that we use to isolate the variable parts of the emission 
lines (e.g., Peterson et al. 2004); our initial impression is that the optical Fell blends do not 
appear in rms spectra and thus there is no compelling evidence for optical Fell variability, 
at least on time scales of several months. 

Thus, the goal of this study is to determine whether or not Fe ll emission varies in NGG 
5548 and, if so, with what amplitude and on what timescale. We also attempt to constrain 
meaningfully the Fe ll line width as an additional test of the common assumption that the 
optical Fell emission arises co-spatially with the Balmer lines (e.g., Phillips 1978). We will 
not specifically re-examine the issue of the variability of the small blue bump since the data 
examined by Maoz et al. (1993) have not been superseded by superior data and we are 
unlikely to be able to improve on their results. We will thus concentrate on variability of the 
optical Fell features, comprised primarily of two blends, one between II7 and H/?, and the 
other covering the spectral ranges 5100-5700 A (often referred to as the 4570 A and 5190, 
5320 A blends, respectively; Osterbrock 1976), and of the UV Fell blends in the vicinity of 
the GlvA1549 emission line. 



2. METHODOLOGY AND DATA ANALYSIS 

2.1. Data 

As noted above, we are focusing on NGG 5548 on account of the large amount of 
variabihty data available. There are two suitable UV data sets. The flrst of these is a 
set of lUE spectra obtained over an 8-month period in 1988-89, sampled once every four 
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days (Clavel et al. 1991). The second set of UV spectra were obtained with the Faint Object 
Spectrograph (FOS) on Hubble Space Telescope (HST) over a 39-day period in 1993, sampled 
once per day (Korista et al. 1995); this second set was supplemented by additional lUE data, 
that we will disregard in this study as they arc far less suitable than the HST data for the 
present purposes. Both UV data sets were recovered from the International AGN Watch 
website. The lUE spectra used were those processed with the lUE New Spectral Image 
Processing System (NEWSIPS; Nichols et al. 1993). The iJST spectra were processed as 
described by Korista et al. (1995). 

In addition to the UV spectra, there arc 1248 optical spectra that were obtained begin- 
ning in 1988 (i.e., concurrently with the lUE spectra described above) for 13 years through 
2001 (Peterson et al. 2002 and references therein). However, only a relatively small subset of 
these data are suitable for studying variability of the optical Fe ll blends; the Fe ll blends are 
broad, low-contrast features that can be measured or modeled usefully only in high-quality 
spectra. Also, as discussed below, we find that Fe ll variability is most reliably measured by 
combining spectra to form difference spectra (i.e., subtraction of a low flux state spectrum 
from a high flux state spectrum) or rms spectra in order to isolate the parts of the spectrum 
that are varying. Since the amplitude of variability is small and since combining spectra 
propagates noise, it is essential to start out with high-quality data. Specifically, we require 
high signal-to-noise ratios, reasonably high spectral resolution in order to detect structure in 
the Fe ll blends, and high-fidehty relative flux calibration over an extended wavelength range 
in order to constrain the underlying continuum level. Moreover, the necessity of combining 
spectra requires that we select a homogeneous subset of the data. Given these requirements, 
the best subset is drawn from those obtained with the Lick Observatory 3-m telescope by 
A.V. Filippenko and collaborators ("Set H" described by Peterson et al. 2002 and references 
therein). We selected from this subset all high-quality spectra that cover the observed wave- 
length range 4000-6000 A. These criteria yielded 73 spectra spanning the entire 13 years of 
the International AGN Watch program in a fairly representative way. These 73 spectra were 
put on a common flux scale by assuming that the narrow-hne [O ill] AA4959, 5007 fluxes are 
constant with F([Ora])A5007 = 5.58 x 10"^^ ergssec^^cm^^ (Peterson et al. 1991). This 
was accomphshed using a version of the van Groningen & Wanders (1992) flux scahng al- 
gorithm. After scaling each spectrum in flux, we verified the consistency of the calibration 
by measuring the continuum flux at 5100 A and the total H/? flux in each spectrum. Figure 
1 shows a comparison of the continuum and Y{j3 fluxes measured here with published fiuxes 
based on all of the International AGN Watch data. Agreement appears to be excellent. 
Figure 1 also underscores the point that the spectra we have selected for study are a suitably 
representative subset of the entire database. 
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2.2. Fe II Templates 

As noted above, the complex atomic structure of Fe+ produces a large number of tran- 
sitions that pervade the UV, optical, and IR spectrum. Doppler broadening within the BLR 
causes the Fe ll emission lines to blend together and with lines of other species, thus making 
it extremely difficult to accurately measure the Fe ll emission- line flux. However, it was first 
noted by Phillips (1977, 1978) that the optical Fell spectra of Seyfert 1 galaxies are very 
similar, differing from one AGN to the next primarily in the amount of Doppler broadening. 
In this context, the Seyfert 1 galaxy I Zw 1 (Sargent 1968; Phillips 1976) is of particular 
interest: the widths of the Fell and Balmer lines in this source are only ~ 1100km s~^, 
among the very narrowest permitted lines observed in Seyfert 1 spectra. Phillips demon- 
strated that by suitably broadening the I Zw 1 spectrum to match the broad H/3 widths of 
other Seyfert Is, the observed Fe ll profiles could be matched reasonably well. This suggested 
that the Fell emission arises co-spatially with the broad H/9 emission. This also suggested 
to Boroson, Persson, & Oke (1985) and Boroson & Green (1992) that the broadened I Zw 1 
spectrum could be used as a template to model the Fe ll contribution to AGN emission-line 
spectra, although of course their interest was more focused on removing Fe ll as an emission- 
line contaminant than in using it to determine the Fell emission flux. Boroson & Green 
used high-quality spectra of I Zw 1 to develop an Fell emission-line template that could be 
used generally to remove Fe II emission from AGN spectra over the range 4250-7000 A. Re- 
cently, a revised and improved optical Fe II template based on I Zw 1 and covering the range 
3535-7534 A has been prepared by Veron-Cetty, Joly, & Veron (2004). It is this template, 
shown in Fig. 2, that we use in analysis of the optical spectra of NGG 5548. Templates 
with a range of intrinsic line widths were generated by folding the entries of Fe ll transitions 
listed in Table A.l of Veron-Getty et al. (their broad system Ll) with Lorentzian^ profiles 
of different FWHM. 

A similar template for the UV spectrum (1250-3090 A), again based on I Zw 1, has 
been constructed from iJST spectra by Vestergaard & Wilkes (2001). We use this template, 
shown in Fig. 3, to attempt to model the UV Fell emission in NGG 5548. Vestergaard 
& Wilkes present several versions of UV iron templates, including separate templates of 
Fell and Felll emission, respectively, each of which have options of including or excluding 
certain multiplets that are characteristically strong in IZwl. Two templates of UV iron 
emission were fit to the NGG 5548 data: (1) a template consisting only of Fell emission 
transitions, excluding Fell UV191 at 1786 A, and (2) a template consisting of both Fell 



^We follow here the authors on the choice of profile shape for the template. However, the actual profile 
shape has little importance in this context, especially once the template is broadened. 
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and Felll emission transitions, excluding Fell UV191 and Felll UV34 as these multiplets 
are particularly strong in I Zw 1 and the spectra of NGC 5548 show little indication of their 
presence. As will be clear from § 2.4, the template that includes Fc ill emission provides a 
poor fit to the data. We note that the Fell emission shortward of the CiV emission line 
tends to be overpredicted by using the template as originally presented by Vestergaard & 
Wilkes. Forster et al. (2001) came to a similar conclusion for a small fraction of the spectra 
in their sample. The bluest part of the iKT spectrum of IZwl (below ~1500A) on which 
the UV template is based was not observed at the same time as the longer-wavelength 
segments and I Zw 1 increased in luminosity between these two observations. Vestergaard & 
Wilkes rescaled the shortest- wavelength segment to match the longer- wavelength part of the 
spectrum to produce a full-range UV template. However, it was recognized then that the 
Fell emission in the template below ~1550A may need to be rescaled for some objects. By 
fitting the NGC 5548 spectra we found that very reasonable fits can be obtained by applying 
a scaling factor of 0.25 to the template Fell emission shortward of 1550 A. 

The general fitting procedure used here is described by Vestergaard & Wilkes (2001; 
section 4.2 in that paper). A power-law continuum is subtracted from the spectrum, which 
makes the Fell emission more visible, before fitting the Fell templates. The power- law 
continuum is fitted to windows that are virtually free of line emission, namely ~4200 A and 
~5650A in the optical and ~1350A and ~2000A in the UV. An iteration over continuum 
setting and Fell template fitting may be necessary to get the best fit. We consider the 
best template fit to be that which leave the residual spectra virtually featureless in regions 
dominated by Fell. Additional details on our fitting of the optical spectrum of NGC 5548 
are given in §2.3.3. 



2.3. Optical Fe II Blends 

We made an initial attempt to fit each of the 73 selected spectra with the optical Fe ll 
template. We found this to be very difficult because of blending of Fe ll with other features 
in the spectra, notably weak narrow lines and small-scale structure due to features in the 
underlying galaxy spectrum. Given the problems we encountered, we had little confidence 
in these results. 



-8- 



2.3.1. Difference Spectra 

We therefore began experimenting with alternative approaches that would mitigate some 
of the problems due to contaminating features. One method is to produce difference spectra, 
i.e., to subtract a low flux-state spectrum from a high flux-state spectrum. This effectively 
removes all the non- variable components from the spectrum, including narrow emission lines 
and host galaxy hght, although sometimes small residuals of these features can be seen, 
usually as a consequence of seeing and aperture effects (e.g., Peterson et al. 1995) or small 
differences in the line-spread function between pairs of observations. A potential liability of 
this approach is that the uncertainties in difference spectra can be fairly large; the differences 
between two spectra are often quite small, but the uncertainties add in quadrature. We 
attempted to mitigate this problem by making a high signal-to-noise ratio low-flux state 
spectrum by co-adding a number of individual faint-state spectra, and then subtracting this 
"standard low-state spectrum" from all of the individual spectra. To form a standard low- 
state spectrum, we identified several spectra obtained during the faintest periods in the 13- 
year monitoring program. The standard low-state spectrum as formed by averaging 3 spectra 
from the fourth year of the AGN Watch program (specifically, the spectra obtained on Julian 
dates 2448733, 2448765, and 2448810) with two spectra obtained during the thirteenth year 
(on Juhan dates 2452030 and 2452045). 

Use of difference spectra greatly improved our ability to model the optical Fe ll emission 
with a template, although the results were still not quite satisfactory. We found that there 
was degeneracy at about the 10% level in fitting the Fell fluxes (i.e., changes in flux of 10% 
made no discernible difference in the overall quality of the flt). Moreover, comparison of the 
inferred Fe II fluxes between observations that are closely spaced in time often showed large 
differences; indeed, if the errors are random, we would infer that the level of precision of 
the Fell fluxes is only ~ 30%, about an order of magnitude more uncertain than our H/3 
fluxes (cf. Peterson et al. 2002). Despite these large uncertainties in individual spectra, we 
nevertheless flnd a clear correlation between the continuum flux at 5100 A and the Fell flux. 
We flnd a Pearson's linear correlation coefficient of 0.73 and a probability of a chance (or 
null) correlation below 0.0001. While we caution about placing too much emphasis on the 
specific numbers, these results emphasize the reality of such a correlation. 

These measurements establish that the optical Fe ll fiuxes do indeed vary in a way that 
is correlated with the continuum. However, we are unable to establish a timescale for the 
response of the optical Fell blends to continuum changes, i.e., to measure a reverberation lag. 
First, the sampling of these 73 spectra is comparatively sparse, as can be seen clearly in Fig. 1. 
The median interval between observations is nearly 60 days. For reference, cross-correlation 
of the 73 optical continuum and H/3 fluxes measured from these spectra and displayed in 
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Fig. 1 yield time delays Tcent = 5.3 ± 11.8 days and Tpeak = 17 ± 4 days, where these 
two quantities are respectively the centroid and peak of the optical continuum/H/3 cross- 
correlation function; both of these measurements are consistent with the result obtained 
using all the AGN Watch data shown in Fig. 1, although the centroid value is formally 
highly uncertain. Second, the very large uncertainties in the optical Fe II fluxes obviate any 
chance of obtaining even a marginal time-lag measurement. Formally, cross-correlation of 
the continuum and Fell light curves yields Tcent ~ 12 ± 27 days and Tpeak = 11 ± 34 days 
with a peak correlation coefficient Vmax — 0.75. We can conclude, however, that the lag for 
the optical Fe ll emission is less than several weeks; this is not an especially stringent upper 
hmit as it exceeds the value of any other rehably measured lag in NGC 5548. 

Our original planned strategy had been to start with the most suitable optical spectra 
and progress to the less suitable, but better time-sampled data. However, we abandoned that 
approach when these initial results made it clear that measuring the Fe ll fluxes accurately 
even in the best data is extremely difficult; measuring the Fell fluxes in less suitable data 
would not improve the situation. 

Inspection of the fits to our 73 difference spectra indicate that the difficulties in accu- 
rately fitting Fell templates to the spectra are largely attributable to noise in the spectra, 
of both Poisson and non-Poisson origin. In other words, wc might hope to obtain more 
accurate flux measurements by combining not only low-state, but high-state spectra as well. 
This proved to be essential for estimating the amplitude of Fell flux variability. In this case, 
we examined the continuum flux measurements for all the AGN Watch data in Fig. 1 and 
computed their mean and standard deviation. We deflned a mean high-state spectrum by 
combining all of the spectra in our selected sample of 73 spectra for which the continuum 
flux is more than one standard deviation above the mean value. These spectra are those 
obtained on Jufian dates 2448243, 2448280, 2449802, 2449830, 2449930, 2449986, 2450336, 
2450361, 2451017, 2451056, 2451077, and 2451221. Simflarly, we deflned a mean low state 
spectrum by combining all of the spectra for which the continuum flux is more than one stan- 
dard deviation below the mean values; the spectra that meet this criterion were obtained 
on Juhan days 2447912, 2448089, 2448102, 2448132, 2448733, 2448746, 2448765, 2448798, 
2448837, 2450664, 2452030, and 2452045. We then computed the difference between these 
two in order to produce a high signal-to-noise ratio difference spectrum, as shown in Figs. 4 
and 5. 
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2.3.2. Mean and RMS Spectra 



Another approach to more clearly isolating the variable part of the spectrum (and thus 
suppressing the constant components that contaminate the optical Fell spectrum) is to 
compute a root-mean-square (rms) spectrum. Using all of the 73 spectra in the sample, we 
formed a mean spectrum, 



where Fi{X) is the ith spectrum of the spectra that comprise the database. Similarly, the 
rms spectrum is defined by 



The mean and rms spectra formed from our 73 optical spectra are shown in Fig. 6. This mean 
spectrum is the only spectrum which we fitted for this discussion that does not strictly contain 
variable emission components only. The mean spectrum was fitted to allow assessment of 
the relative Fe ll variability ampfitude. 

2.3.3. Identifying the Best Fitting Method of the Optical Iron Emission in NGC 5548 

The 4450-4650 A and 5050-5600 A regions of AGN and quasars are commonly dom- 
inated by Fell emission which can be reasonably well-fitted by an optical Fell template 
generated from the I Zw 1 emission spectrum (Boroson & Green 1992; Veron-Getty et al. 
2004). However, we were unable to fit these regions in the NGG5548 optical spectra with 
either the Veron-Getty et al. or the Boroson & Green templates in their original form, i.e., 
with relative line intensities as observed in I Zw 1 . The top panel of Fig. 7 illustrates the 
result of accounting for all the fiux at ~ 4500 A above a power-law continuum with the 
Veron-Getty et al. Fell template. Glearly, this strongly overpredicts the Fell fiux in the 
5050-5700 A range. Fitting the template to match the 5050 - 5600 A Fe ll features instead 
leaves strong residuals at ~ 4500 A. 

We examined three possible solutions to this problem. First, the Fe II features may not 
have the relative line intensity as observed for I Zw 1. Second, the emission at ~4500 A may 
be due to extended wings of some combination of H7, He II A4686, and H/3. Third, additional 
line emission, previously unaccounted for, may be present. We address each in turn in the 
following. 




(1) 



i=l 
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1. We did indeed find that the spectra could be fitted well if the 5050-5600 A Fe ll features 
in the template arc rescaled. The best fits to the mean, difference, and rms spectra 
of such modified Veron-Cetty et al. (or Boroson & Green) templates yields scalings 
relative to (fixed) ~4500A Fell strength of ~0.4, 0.25, 0.2, respectively. However, 
such low intensity ratios are not commonly seen in AGN spectra; relative intensities of 
0.7 - 0.8 seem more typical. Also, it is apparently difficult to obtain very low intensity 
ratios of these fines in standard photoionization equilibrium modeling (M. Joly 2004, 
private communication). This solution is hence not particularly attractive. 

2. To test if extended wings of the H7, Hell, and H/5 lines can account for the ~ 4500 A 
residual emission seen in the lower panel of Fig. 7 we fitted Gaussian functions to H7 
and He ll A4686, while a combination of Gaussian and Lorentzian functions was neces- 
sary to reproduce the (variable) H/3 profile. These modeled lines were then subtracted 
from the data. The results are shown in Fig. 8. The upper panel shows the best-fit 
width to the Hell line of 10,000km s-^FWHM). A strong residual feature at 4500 A 
shows that we still cannot account for all the fiux in this region. A second attempt 
with a broader (12,000km s"^ FWHM) Hell line, as shown in the lower panel of Fig. 
8, results in very little improvement. 

3. An alternative explanation is that additional previously unidentified line emission is 
present in the NGC 5548 spectrum. Although emission lines other than Fe ll are rarely 
seen at ~ 4500 A in quasars (e.g., Boroson & Green 1992), NGC 5548 may have a strong 
contribution of He I A4471, given the presence of strong He I A5876 observed at the edge 
of the spectra (Figs. 4, 5, and 6). In fact, the presence of HelA5876 also implies the 
existence of the He I singlet fines at 4922 A and 5016 A. Common intensity ratios to 
He I A 5876 are 0.3, 0.1, and 0.1 for the 4471 A, 4922 A, and 5016 A lines, respectively, 
although the 4471 A intensity ratio may be as low as 0.1 (K. Korista 2004, private 
communication). Both panels of Fig. 8 illustrate that the excess emission at ~ 4500 A 
remaining when subtracting other known sources of line emission is consistent with the 
presence of He I A4471 (thin solid curve) at an intensity of 0.3 times that of He I A5876. 

These arguments collectively justify our inclusion in our general fitting of the non- 
apparent He I AA4471, 4922, 5016 fines with the above-mentioned intensity ratios as found in 
photoionization equilibrium modeling. Figure 5 illustrates the method we thereafter adopted 
for fitting the Fell emission: the expected He I features, based on the observed 5876 A 
line, were subtracted before fitting the Fell emission. The adopted template fitting of the 
difference, mean, and rms spectra are shown in Figs. 4 and 6. 

We note that both the broad He I A5876 and He ll A4686 features are observed shifted to 
the blue by ~1000km s~^ to 1500 km s~^; this is apparent even in the observed spectra. It 
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is particularly clear in the mean spectrum (Fig. 6) where the narrow He I A5876 component 
is observed at 5876 A yet the peak of the broad component is blueshifted; we reproduced 
this shift when modeling and subtracting the expected He I lines in the mean spectrum. 
Interestingly, there is no clear indication that the (broad) HellA1640 line is shifted by such 
a large amount. 

Upon subtraction of the best-fit Fell templates, a few small residuals may still remain, 
especially in the mean spectrum (Fig. 6) where some non- variable possibly narrower Fell 
transitions may be present; we did not attempt to include any transitions of the narrower 
Fell system identified by Veron-Cetty et al. (2004), on the assumption that these would 
not appear in our difference or rms spectra, which appears justified by Figs. 4-6. It 
is noteworthy that many of the smaller-scale features seen atop the 4450-4750 A emission 
blend in the mean spectrum are consistent with the expected positions of Fe ll m38. 

2.3.4. Results 

Our analysis of the optical Fell data on NGC 5548 leads us to the conclusion that the 
Fell emission is indeed variable in flux and that these variations are positively correlated 
with continuum variations. We are unable, however, to make any meaningful determination 
of any time delay between the continuum and Fe ll emission blend variations. 

In an attempt to assess the amplitude of variability of the optical Fe ll emission, we 
integrated the Fe ll flux in the best-flt template over the range 4250-5710 A and divided it 
by the H/3 broad-line flux as measured in the mean, difference, and rms spectra, respectively. 
The results are shown in Table 1. Comparison of the Fe ll/H/3 flux ratio in the variable spectra 
with the mean spectra shows that this ratio is only 25-50% smaller in the variable spectra. 
This implies that the optical Fell emission varies with an amplitude of 50-75% that of the 
broad H/? component. The reason that Fell emission is not always obvious in rms spectra is 
simply because the variations are generally rather small, and the low-contrast Fell features 
are weak relative to H/3 and do not show up very well even in high-quality difference or rms 
spectra. 

The best fit we achieve is with an optical Fell template broadened to a width similar 
to that of the variable H/3 line; the HjS line width in the rms emission-line spectrum shown 
in the lower panel of Fig. 6 is about 6250 km s~^. Interestingly, this template gives only a 
slightly improved fit to the 5050 - 5600 A Fell features relative to templates of somewhat 
broader widths of ~10,000km s~^; the differences in the fit and residuals are practically 
insignificant. 
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It is also worth noting, however, that for broadening larger than ~ 9000 km s ^, the 
templates become essentially degenerate, i.e., it is hard to distinguish among them. 

Given that reasonable fits are obtainable with a range of template line widths, we carried 
out additional tests to try to determine the maximum and minimum line widths for which 
the template fitting process gives unacceptable results. We find a range of template fine 

widths that yield "acceptable" fits, and a smaller range that are "preferred," for which even 
better fits are obtained. The degeneracy is mainly due to the flux density uncertainties in 
the difference spectra (section 2.3.1) and the relatively minute differences in the detailed 
shapes of Fell templates that differ by 250-750 km s~^ in line widths, especially for large 
line widths (above about 6000 km s^^). For the optical difference and rms spectra, which are 
not contaminated by non-variable emission, we find that the template fits are unacceptable 
at line widths less than 4500 km s~^; "preferred" fits were obtained for widths above 5750 km 
s~^, and our "best fit" value of 6250 km s~^ is a subjective assessment that this provides 
the best overall fit to the spectrum. Establishing an upper limit on the broadening is more 
difficult because of the stronger degeneracy of the templates with line widths above about 
9000 km s~^ We are unable to place any meaningful constraints on the upper limit to the 
Fe II line widths. 

Subtraction of the optical Fe ll template leaves He ll A4686 in the residual spectrum. The 
width of this line in the residual spectrum is very broad, ~ 12, 000 km s~^, and the width is 
insensitive to which Fell template is used, i.e., the difference between using the 6250km s~^ 
and 10,000 km s^^ templates changes the HeIlA4686 width by less than about 200 km s~^. 
This HellA4686 line width is very similar to the width of HellA1640 measured in the Fell 
subtracted UV difference and rms spectra, discussed below, and to the HellA1640 width 
obtained during both UV monitoring campaigns, and significantly larger than the width 
of He II A4686 measured in the rms spectrum, with no attempt to account for the possible 
effects of Fell, during the first year of the NGC 5548 monitoring campaign (Peterson et al. 
2004). It is not clear, however, whether the HellA4686 width varies with time, as does, 
for example, the width of H/9, or if the base of the HeIlA4686 line gets broader when the 
Fell is removed. Certainly, the responses of HeIlA1640 and HeIlA4686 are expected to be 
identical (Bottorff et al. 2002), and the lags are, in fact, consistent within their uncertainties 
(Peterson et al. 2004). Indeed, there seems to be no strong evidence that either the lags or 
line widths of HellA1640 and HellA4686 are significantly different from one another. The 
concerns of Bottorff et al. about the differences in the response of these two fines may be 
not as significant as previously thought. 
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2.4. UV Fell Blends 

2.4-1- UV Difference Spectra 

In order to determine whether or not the UV Fe ll emission blends vary with the contin- 
uum, we carried out an analysis similar to that described above for the optical emission lines, 
though in this case we describe only our investigation of difference spectra. For both UV data 
sets described in §2.1, we computed mean high-state and mean low-state spectra, again using 
the criteria that these spectra have continuum fluxes that are one standard deviation above 
and below the mean flux, respectively. The HST high state spectrum included those obtained 
on Juhan dates 2449121, 2449122, 2449123, and 2449124, and those that comprise the low- 
state spectrum are those obtained on 2449097, 2449098, and 2449100. These spectra and 
their difference are shown in Figs. 9 and 10. For the lUE data, the high-state spectrum was 
constructed from the spectra obtained on Juhan dates 2447538, 2447543, 2447546, 2447621, 
2447625, 2447629, 2447633, 2447634, 2447637, 2447638, 2447641, 2447645, and 2447649, and 
the low-state spectrum is based on the spectra obtained on Julian dates 2447577, 2447581, 
2447688, 2447692, 2447729, 2447733, 2447737, 2447741, and 2447745. These spectra and 
their difference are shown in Figs. 11 and 12. 



2-4-2- Results 

We fit the difference spectra with two types of template, (1) a pure Fell template, and 
(2) a template that includes both Fell and Felll. Both templates are from Vestergaard 
& Wilkes (2001) and described in § 2.2. While Fell is clearly present in the difference 
spectrum, it is not clear that Felll is present. If Felll is present, it is very weak. We flnd 
that inclusion of Felll is not necessary to model the spectrum, and the spectrum is better flt 
by the pure Fell template. One possible exception is the bump at ~2200A, which is most 
likely Felll (Vestergaard & Wilkes 2001). However, the results of this study are not sensitive 
to the strength of this feature and we will not consider it any further. Compared to the 
HST spectrum, the lUE spectrum is harder to fit because it is noisier and covers a shorter 
wavelength range. On the other hand, if we constrain the slope of the /C/E difference spectrum 
by using the slope measured in the HST difference spectrum, we obtain a satisfactory fit. 
We find that a better fit is obtained with the 10,000 km s^^ template than for the 6250 km 
s^^ template (Figs. 10 and 12). In particular, we find unacceptable fits using templates with 
broadenings less than 8250km s~^. Given the same issues pertaining to the optical spectra, 
we are unable to place any upper limits on template broadenings which produce acceptable 
fits. 
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We also compare the strengths of the UV Fe ll blends (integrated across the two intervals 
1430-1 500 A and 1685-1850 A) and H/3 in the variable part of the spectrum with what is 
observed in the mean spectrum, as we did with the optical Fell emission above^. The results 
of this comparison are given in Table 2. We see here that the UV FeIl/H/5 flux ratio is 
significantly larger in the variable part of the spectrum than in the mean spectrum (we 
note in passing that analysis of the UV mean and rms spectra leads to the same conclusion) . 
Unlike the optical Fe ll emission, the UV Fe ll emission seems to vary with a larger amphtude 
than H/3. 

3. DISCUSSION 

The principal result of this study is that the optical and UV blends of Fe II emission in 
the spectrum of NGC 5548 do indeed vary with time. These variations are correlated with 
the continuum variations. The optical Fell blends vary with an amphtude 50-75% that of 
H/3, and the UV Fell blends seem to vary with an even larger amplitude than H/3. The 
variations of these low-contrast features are most clearly seen in difference or rms spectra, 
which remove non- variable structure in the spectra. However, the variable Fell features are 
not always obvious in difference or rms spectra simply because these features are often lost 
in the noise. Previous claims of the lack of variability of Fell seem to be consistent with 
this conclusion, given the sensitivity of the observations. Unfortunately, on account of the 
weakness of these features in difference or rms spectra, we are unable to determine the time 
scale for response of the Fe ll emission lines to continuum variations. 

We note in passing that very broad Fell features produce a "pseudo-continuum" (e.g.. 
Wills et al. 1985) that is not readily identified as line emission. For example. Figs. 9 and 
11 (upper panels) show that the 1700-1900 A region in NGC 5548 is relatively flat, giving 
the appearance that no (variable) emission lines are present, even though this region is at a 
somewhat higher level than the continuum at ~1450A and 2025 A. However, the presence of 
Fe II emission in this spectral range is clear when fitting the spectrum with an Fe ll template. 

Interestingly, we find that the UV Fell features in the difference or rms spectra are 
best fitted with templates that have a large velocity width, ~10,000km s~^; the optical Fell 
features are consistent with both such a broad template and one of the variable H/3 line 
width. We have no ready explanation for this result, although we might speculate that it 
may be related to the issues responsible for the difficulty of accounting for the strength of 



^The H/3 fluxes were measured in mean or diff'erence spectra constructed from optical spectra that are 
contemporaneous with the UV spectra used. 
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the Fell (Baldwin et al. 2004), including the unusual smoothness of the Fell emission; either 
microturbulcncc or velocity shear seems necessary to account for the discrepancy between 
theory and observation. 

The issue of whether or not the optical and UV Fe ll is emitted from the same region 
is important for theoretical modeling of their relative emission strengths (Baldwin et al. 
2004). While the errors on the Fell emission- line response time are very large, it is possi- 
ble to constrain the size of the Fe Il-emitting regions indirectly through the line-broadening 
parameters under the assumption that Fell follows the same virial relationship seen in ev- 
ery other wcll-charactcrizcd variable line in NGC 5548 (Peterson & Wandel 1999). If we 
assume that the virial product CTcentFWHM^/G is the same for all lines in NGC 5548, using 
measurements from Peterson et al. (2004), we can conclude that the optical Fen emission, 
with FWHM > 4500km s~^, arises at r < 24 hght days, and the UV Fell emission, with 
FWHM > 8500 km s~^, arises at r < 7 light days. This hmit on the size of the optical 
Fell emission gas does not violate the weak constraint on the time delay of several weeks 
based on cross- correlation of the continuum and optical Fc ll emission light curves (section 
2.3.1). Given the assumptions, this is a weak conclusion, but it is worth noting that dif- 
fering line widths for the UV and optical Fell emission suggests that the UV and optical 
Fe II might be emitted predominantly from different regions within the BLR. If indeed the 
BLR gas density increases steadily and smoothly toward the center of the BLR, the distance 
limits are broadly consistent with expectations from photoionization models (e.g.. Fig. 11 of 
Baldwin et al. 2004), namely that UV Fell is most efficiently emitted from higher-density 
gas (presumably located closer to the central source), while optical Fell is more efficiently 
emitted from lower-density gas farther out. 

Veron-Cetty, Joly, & Veron (2004) find that the I Zw 1 Fe ll emission spectrum is com- 
prised of two components of different widths, ~1100km s~^ and ~300km s~^, the latter 
mostly being forbidden iron transitions. This suggests that some of the Fe ll emission origi- 
nates in the narrow- line region. The fact that no variations are detected in the narrow Fe ll 
features by virtue of their absence in the residual spectra of the difference and rms spectra 
is consistent with this explanation; blends of narrow Fc II lines may be present in the mean 
optical spectrum, especially longward of [Olll]. We cannot identify narrow Fell in any of 
the UV spectra we have examined. 

Our conclusion that the amplitude of optical Fe II variability is somewhat smaller than 
that of H/3, but is often of too low contrast for reliable detection, seems consistent with results 
for other sources (KoUatschny et al. 1981; KoUatschny & Fricke 1985; Giannuzzo & Stirpe 
1996; Doroshenko et al. 1999). Apparent exceptions are Fairall 9 and NGC 7603. KoUatschny 
et al. (2000) found the optical Fe ll to vary as strongly as H/3 for NGC 7603. In the case of 
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Fairall 9, KoUatschny & Pricke find that both H/3 and the optical Fell emission vary with 
large amplitudes, by factors of 5.6 and 2.6, respectively. For more direct comparison with 
this result, we restricted our Fell measurement window to the 4570A, 4928 A, and 5018 A 
multiplets, the approximate region used by KoUatschny & Fricke, and found similar values 
to those we quote in Table 1. We conclude that the relative amplitudes of Fell and H/3 flux 
variability amplitudes can differ significantly among AGNs. 

As current photoionization models fall short in reproducing the observed AGN Fell 
emission, Collin & Joly (2000) suggest that Fell in AGNs is due entirely to shock (i.e., 
collisional) excitation. They draw parallels with strong Fell emitting stellar sources which 
have no X-ray or ionizing emission, but do have strong variability and strong outflows. 
Baldwin et al. (2004) do find that collisional excitation can explain the observed Fe ll, though 
they are reluctant to consider introduction of another component, which produces no other 
significant emission, into an already complex picture of the BLR. Moreover, it is not clear 
how such a model could reproduce the observed correlation between the optical Fe ll and the 
continuum fluxes that we observe. 

The template flts show that the IZwl Fell templates are generally reasonably good 
matches to the optical and UV Fell emission from NGC5548. The exceptions are, as noted 
earlier, the apparent absence of the Fe III transitions over the range 1825-2100 A and the poor 
match of the original templates at 2025 A where the Fe ll emission appears to be missing. 

4. CONCLUSIONS 

Our chief conclusions of our study of the Fe ll emission in NGC 5548 can be summarized 
as follows: 

1. The optical and UV Fen emission blends do indeed vary in response to variations 
in the continuum emission, i.e., the Fell and continuum flux variations are strongly 
correlated. Detection of these variations is difficult because these lines are very broad 
and of low contrast. 

2. The optical Fell varies in strength with an amphtude of 50% - 75% of that of H/3. 
The UV Fe ll varies much more strongly. 

3. These results are consistent with most earlier studies claiming no variability owing to 
the difficulty of detecting small-scale variations in low-contrast features. 

4. Unfortunately, the weakness of the Fell features preclude a meaningful variability time 
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scale of this emission to be obtained. The flux variations of the optical Fe ll blends lag 
behind those in the continuum by no more than several weeks. 

5. For reasons that are unclear, we flnd somewhat improved flts using UV Fell templates 
with large widths of order 10,000 km s"^, as opposed to templates with widths similar 
to that of H/3. 

6. Fell removal improves measurement of other emission lines. The HellA4686 feature in 
the iron-subtracted spectrum displays a line width consistent with that of Hell A1640, 
as expected from photoionization models. Failure to account for the effect of Fell 
emission seems to yield a He ll A4686 line width that is too small. 

7. Because we cannot meaningfully constrain either the line width or the reverberation 
response timescale for Fell, we are unable to determine conclusively whether or not 
the Fe ii emission originates in the same region as H/5 or whether the UV and optical 
Fe II is spatially co-emitted; this issue has importance for the theoretical modehng of 
the relative UV and optical Fell emission strengths (e.g., Baldwin et al. 2004). 

However the line widths below which the optical and UV templates produce unac- 
ceptable flts to the observed Fell emission do suggest that the UV and optical Fell 
may originate in different regions, although there may be some overlap of these regions 
given our inability to place meaningful upper limits on the Fe ll line widths. 
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Fig. 1. — Optical continuum (upper panel) and H/3 (lower panel) light curves for NGC 5548 
for the period 1988 to 2001. The vertical lines show the ±1(7 error ranges for the hght 
curves based on the International AGN Watch database (Peterson et al. 2002 and references 
therein), and the open and filled circles show the measurements based on the 73 rescaled 
spectra used in this study. The filled circles mark the epochs of the data used to construct the 
low and high state spectra shown in Figure 4 and described in section 2.3.1. Continuum fluxes 
are in units of 10^^^ ergs s^^ cm^^ A^^, and H/? fluxes are in units of 10^^'^ ergs s^^ cm^^. The 
distribution of open and filled circles shows that this particular subset of spectra represents 
the entire database very well. 



-23- 



I 

CM 
I 



S 0.4 - 



o 

I 

u 
0) 



I 

o 



0.2 




0.4 - 





H — I — I — — I — I — ^ 



H — I — h 



\ 1 1 1 1 H 





4000 5000 6000 7000 

Rest Wavelength (A) 



Fig. 2. — The upper panel shows the optical Fell template constructed by Veron-Cetty, Joly, 
& Veron (2004), reproduced as described in §2.2; the intrinsic line width is 1000 km s~^ 
(FWHM). The lower panel shows this template with an intrinsic Doppler width of 6250 km 
s-i (FWHM). 
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Fig. 3. — The upper panel shows the unbroadened ultraviolet Fell template constructed 
by Vestergaard & Wilkes (2001). The fluxes shortward of 1550 A are scaled by a factor 
0.25 compared to the original template of Vestergaard & Wilkes (see §2.2), which yields a 
better fit in general to the data analyzed here. Note that no Fe ill emission is present in this 
template. The lower panel shows this template artificially broadened to a Doppler width of 
6250 km s"^ (FWHM). 
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Fig. 4. — The upper panel shows the high state and low state optical spectra of NGC 5548 in 
the rest frame. The lower panel shows the difference spectrum formed by subtraction of the 
low state spectrum from the high state spectrum (upper spectrum in the upper panel) along 
with our preferred final fit to the spectrum, as described in §2.3.3.. Note that the residuals 
in the [O ill] AA4959, 5007 hues arise from small differences in the hne-spread function; the 
total line flux in the residuals is nevertheless zero. The smooth solid thick hue shows the 
broadened (FWHM = 6250km s~^) and scaled Fell template, the dashed line shows a 
fit to the underlying continuum, and the dash-dot line shows a Gaussian fit (~12,000km 
FWHM) to the residual HeIlA4686 profile. The smooth solid thin line represents the 
He I AA 4471, 4922, 5016, 5876 lines; for visibility, these lines are shown on top of the Fell 
template fit. The spectrum after subtraction of the Fell template, the He I lines, and the 
He II fit is also shown (dashed line) . For visibility a section of this residual spectrum and the 
underlying continuum is shown shifted to a lower flux level. 
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Fig. 5. — Illustration of our method for fitting the optical spectra of NGC 5548. An enlarge- 
ment of the difference spectrum from the lower panel of Fig. 4 with the continuum subtracted 
is shown in the upper panel. The He I AA 4471, 4922, 5016, 5876 fines (dashed) and the spec- 
trum remaining after their subtraction (dot-dashed) are also shown. The AA4471, 4922, 
5016 profiles are those expected based on the observed A 5876 line (see text). In the lower 
panel the smooth solid curve shows the 6250km s~^ broadened and scaled Fell template su- 
perposed on the continuum subtracted and He I subtracted spectrum from the upper panel. 
The lower of the two solid histograms is the residual from subtracting this iron template fit. 
A simple (10,500 km s^^) Gaussian fit to the HellA4686 profile (dot-dashed) and residuals 
with the Fell template and Hell fit subtraction is also shown. Weak residuals remain at 
4543 A. 
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Fig. 6. — The mean (upper panel) and rms (lower panel) optical spectra of NGC 5548 in 
the rest frame. In both cases, the smooth solid thick line shows the broadened (6250 km 
FWHM) and scaled Fe ll template, the dashed hne shows a fit to the underlying continuum, 
and the dash-dot line shows a Gaussian fit to the residual He ll A4686 profile. As in Fig. 4, 
the thin solid fines atop the Fell template fit show the four He I AA 4471, 4922, 5016, 5876 
fines. The spectrum after subtraction of the Fe ll template and the He ll fit is also shown. 
In the upper panel a section of this residual is shown shifted downward for visibility. The 
He II profile in the mean spectrum can be fit with two Gaussian profiles of widths ~800 km 
s~^ and 12,000km s~^, respectively. The width of the rms Hell profile fit is 12,250km s~^. 
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Fig. 7. — Sample fit with the Vcron-Cetty ct ah (2004) template to the difference spectrum 
shown in Fig. 5. In the upper panel, the template is fitted to match the 4500 A region only. 
Notice the strong overprediction of the visible Fe ll features between 5050 A and 5550 A. The 
lower panel shows the template fit and residuals obtained if the template is fit to the fiux 
level of the ~5300A Fell transitions only. This template poorly matches the spectrum in 
the 4500 A region. The Fell template has been broadened to 6250km s . 
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Fig. 8. — The continuum and Fe Il-subtracted difference spectrum from the lower panel of 
Fig. 7 is shown here with the residuals (thick line) obtained when subtracting reasonable 
reproductions of the non-Fe profiles (dot-dashed lines) . The profile fits to H/3 and H7 were 
generated to reproduce especially the blue and red sides of the core and the broad wings as 
allowed by the data; the details of the fits to the hne peaks can be ignored and the fit to 
the [Olll] residuals is suggestive only. The HenA4686 profile is shown with a 10,000 km s~^ 
Gaussian profile (best fit) in the upper panel, while the lower panel shows the widest Hell 
profile of 12,000 km s^^ allowed by the data. Neither fit can adequately account for the excess 
emission at 4500 A. However, reasonable consistency is found with the expected strength of 
HclA4471 based on the observed HelA5876 feature. Here HelA4471 is represented by a 
Gaussian fit to HelA5876, shifted to 447lA (with fixed FWHM) and rescaled by a factor 
0.3, a commonly established strength of HelA4471. This justifies our inclusion of the He I 
fines in our fitting of this region (e.g.. Figs. 4, 5, and 6). 
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Fig. 9. — The upper panel shows the high state and low state ultraviolet spectra of NGC 
5548 in the rest frame as observed with HST in 1993. The lower panel shows the difference 
spectrum formed by subtraction of the low state spectrum from the high state spectrum. The 
smooth solid line shows the broadened (10,000km s^^ FWHM) and scaled Fell template, 
and the dashed line shows a fit to the underlying continuum. 
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Fig. 10. — An enlargement of the difference spectrum from the lower panel of Fig. 9 with the 
continuum subtracted. Again, the solid line shows the broadened and scaled Fe ll template, 
which has a width of 6250 km s^^ in the upper panel and of 10,000 km s~^ in the lower panel. 
Note the poor fit of the 6250 km s^^ iron template to the 1700 A region. The template 
overpredicts the Fell emission at 2025 A, possibly due to strong multiplets present in IZwl 
but not in NGC 5548. The emission on either side of the expected Fe ll bump at 2025 A is 
almost certainly Felll (see Vestergaard & Wilkes 2001). 
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Fig. 11. — The upper panel shows the high state and low state ultraviolet spectra of NGC 
5548 in the rest frame as observed with lUE in 1988-89. The lower panel shows the difference 
spectrum formed by subtraction of the low state spectrum from the high state spectrum. The 
smooth solid line shows the broadened (10,000km s~^ FWHM) and scaled Fell template, 
the dotted hue represents the Fe ll subtracted spectrum, and the dashed hne shows a fit to 
the underlying continuum. 
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Fig. 12. — An enlargement of the difference spectrum from the lower panel of Fig. 11 with the 
continuum subtracted. Again, the solid hne shows the broadened and scaled Fe ll template, 
which has a width of 6250 km s~^ in the upper panel and of 10,000 km s~^ in the lower panel. 
As in Fig. 10, the template broadened to the rms width of broad H/3 (6250 km s~^) yields a 
less satisfactory fit than that of the 10,000 km s~^ template, especially at about 1700 A. The 
Fe II subtracted spectra are shown as the (lower) thin-lined histograms. 
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Table 1. Optical Fell/H/? Flux Ratio 



Template Width 


Mean 


Difference 


RMS 


(km s~^) 


Spectrum 


Spectrum 


Spectrum 


(1) 


(2) 


(3) 


(4) 


6250^ 


0.66 


0.48 


0.30 



^Template fits with 10,000 km s ^ broadenings yield very 
similar results. 



Table 2. Ultraviolet Fell/H/3 Flux Ratio 



Template Width 


HST Mean 


HST Difference 


WE Mean 


lUE Difference 


(km s~"^) 


Spectrum 


Spectrum 


Spectrum 


Spectrum 


(1) 


(2) 


(3) 


(4) 


(5) 


6250 


0.45 


3.12 


0.56 


4.98 


10000 


0.51 


3.12 


0.67 


1.66 



